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Abstract. A 850 /im map of the interacting spiral galaxy M 5 1 shows well-defined spiral arms, closely resembling the struc- 
tures seen in CO and HI emission. However, most of the 850 /im emission originates in an underlying exponential disk, a 
component that has not been observed before in a face-on galaxy at these wavelengths. The scale-length of this disk is 5.45 
kpc, which is somewhat larger than the scale-length of the stellar disk, but somewhat smaller than that of atomic hydrogen. Its 
profile can not be explained solely by a radial disk temperature gradient but requires the underlying dust to have an exponential 
distribution as well. This reinforces the view that the submm emission from spiral galaxy disks traces total hydrogen column 
density, i.e. the sum of H2 and Hi . A canonical gas-to-dust ratio of 100±26 is obtained for Kgso ~ 1-2 g^^cm^, where ksm) 
is the dust opacity at 850 nm. 



1. Introduction 

After the introduction of the SCUBA bolometer camera at 
the JCMT a number of nearby galaxies have been imaged in 
the subm m range to inves t igate the cold dust component of 
the ISM ( Alton et al."2002^ Israel et al. ' 1999: 'Stev ens & Geail 



[2000). The spatial distribution of dust and its relation to other 
components of the ISM is best studied in face-on systems, 
but at submm wavelengths surface brightnesses are typically 
low. This particularly hampers the search for diffuse emis- 
sion from inter-arm regions. Moreover, most of these nearby 
galaxies are larger than the field-of-view of SCUBA requiring 
scan-mapping which, as explained below, causes features in the 
background at the level of a few milli-Jansky frustrating further 
attempts to determine the level of any weak extended emission. 
The example of NGC 6946 is illustrativ e. Little diffuse em is- 
sion in the inter-arm region was seen bv lAlton et al] ( I2OO2I) in 
the 850 /im map. Consequently, they state that "relatively little 
is known about how the dust is distributed with respect to the 
spiral arms and, in particular, whether inter-arm grain material 
is prevalent or not". 

At resolutions much lower than those provided by the 
JCMT, the cold dust in galaxies has been mapped with 
ISO in the far -infrared ( Haas et al. 'l998'; Popescu et al. '2002'; 
Hippelein et aL.2003: Tuffs & Gabriel 2003: Popescu & Tuffs. 
2003b . These observations confirm the ubiquitous presenc e of 
cold dust below 20 K first observed bv lChini et al.l (fl986') and 
indicate that the cold dust component is smoothly distributed 
over the disk and heated by the diffuse interstellar radiation 



llHiDt?elein et al]l2003l: IXu et aljfl994l) . lAIton et alJ tOQli re- 
ported evidence for diffuse inter-arm dust in SCUBA observa- 
tion s of NGC 7331 . For a recent discussion of 'dusty disks' see 
aIso lBianchil ( l2004 . 

M 51 is a well studied, face-on interacting spiral galaxy 
at an assumed distance of 9.7 Mpc. It is well known for its 
strong CO emission, but contrary to NGC 6946 its density- 
wave and spiral arms are well-organized and prominent, show- 
ing a hi gh contrast i n the radio continuum (see e.g. Tilanus 
& AIIen ll99ltll989l) . The distribution of cold dust in the nu- 
clear region of M 5 1 was previously mapped at 1.2 mm by 
lGnelinetalJlll994 . who noted a close coiTelation with CO 
line emission. Our cuiTent 850 /im observations (Fig.0 show 
that its dust emission is exceptionally strong. The spiral arms 
are clearly distinguished, but there is also a substantial amount 
of emission from an extended diffuse disk. To our knowledge, 
such a submillimeter disk has not been clearly identified in any 
other galaxy, but based on the ISO observations discussed pre- 
viously, it may be a common component of spiral galaxies. In 
this paper, we will explore the nature and properties of this disk. 
The 850 fim spiral structure of M 5 1 will be discussed in a fu- 
ture paper. 

2. Observations and data reduction 

M 51 was observed at 850 fim and 450 fim in the spring of 
1998 and 1999 at the JCMT^ using the Submillimeter Common 
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User Bolometer Array (SCUBA: [Holland et alJ ( ll999l) V In or- 
der to map a 13.5' x 13.5' region around M 51, the camera 
was used in scan-map mode during which it is scanned at a 
rate of 3" per second across the field, while the secondary was 
chopped with a frequency of 7.8 Hz in right ascension or in 
declination to cancel atmospheric signal variations as well as 
bolometer DC fluctuations. This uses the revised version of 
the E merson-Klein-Haslam algorithm proposed by Emerson 
ifl995h . The scan-angle across the field, (15. 5 ±60°), was cho- 
sen such that the resulting image is fully sampled. The total 
integration time was about 20 hours spread over 6 nights. 

To restore the brightness distribution of the source, the chop 
must be deconvolved from the observed map using a division in 
the Fourier plane ( Jenness et al. 2000). The Fourier-transform 
(FT) of the chop is a sine wave with zeroes at the origin and at 
harmonics of the inverse chop throw. The initial observations 
used chop-throws of 20, 30, and 65 arcsecs. To minimize re- 
gions of low weight (close to zeroes) in the Fourier plane we 
later added throws of 44 and 68 arcsecs. Each individual obser- 
vation was flat-fielded, corrected for atmospheric opacity, and 
despiked. For each of the 10 different chop-configurations the 
individual maps were coadded. Next, the resulting 10 maps 
were Fourier transformed, weighted, and coadded. The coad- 
ded image was transformed back to yield the finished image. 
This is the standard reduction for SCUBA scan-map observa- 
tions as implemented in the SCUBA User Reduction Facility 
software (SURF: Jenness & Lightfoot ll998l) . 

The data were flux-calibrated using Mars and the sec- 
ondary standards HL Tau and CRL 618. In the final image, 
the beam-size at 850 /im is about 15" (FWHM). The caH- 
brated image is shown in Fig. Q and has an overall rms of 
^9 mJy/beam. The inconspicuous source to the north is the nu- 
cleus of NGC 5 1 95 . The familiar spiral pattern is clearly visible 
and can be traced over a large fraction of the disk. This is quite 
unlike the submm continuum observations of other face-on spi- 
rals, which in general lack the sensitivity to trace the arms out- 
side the inner region. In detail, the morphology of the arms re- 
sembles closely the one seen in CO emission ( Aalto et al. 1999; 
IWielebinskietalJl999l) . and we will address the contamination 
of the continuum emission by in-band line emission in section 
3. 

Remarkably, the 850 fim emission from M 51, showing the 
distribution of the cool dust, is dominated by a diffuse exponen- 
tial disk, a component that has not been directly imaged before 
at sub-mm wavelengths and with a high resolution in a face-on 
spiral. The 450 /im observations also show spiral structure and 
an extended disk. In part because of the difficulty to derive an 
accurate flux calibration and the poorer image quality, we will 
not further discuss the 450 fim data in this paper 

2. 1. Sky noise and background removal 

Analyzing extended structures in SCUBA scan maps is a del- 
icate process. While there is no doubt about the reality of the 
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Fig. 1. SCUBA calibrated and background corrected 850 /im 
grey-scale image of M 51 with 850 fim contour superimposed. 
The contour levels are 13.5 (=1.6 a), 27, 40.5, 58.5, 72, 90 and 
117 mJy/beam. 13.5 mJy/beam corresponds to 2.6 MJy sr~^. 
This image clearly shows the extended 850 /im disk. The small 
source to the north is the nucleus of NGC 5195. 

exponential disk, careful processing is required for realiable pa- 
rameter extraction. There are a number of reasons for this. 

Firstly, while regions of low weight have been minimized 
through the use of many different chop configurations, the zero 
near the origin of the Fourier plane can not be suppressed. This 
results in a poor determination of the largest-scale features in 
the map, including the total flux-density. The effect is similar 
to that caused by missing short-spacings in radio interferometry 
observations and it has a similar effect on the background. 

Secondly, chopping removes noise due to sky fluctuations 
with a frequency comparable to the chop-frequency. It cannot 
remove the effect of a sky signal which is steadily increasing 
or decreasing while scanning from one side of the field to the 
other, i.e. a slowly and monotonically varying component of 
the sky. In such situations the 'on' minus 'off will always 
leave a small residual sky signal. We have used the standard 
SURF option, scanj-lb, to subtract a linear baseline from each 
bolometer-scan fitted to the extreme 2 arcmin on each side 
of the scan. While this significantly improves the flatness of 
the background, it will leave any ripples that may be due to 
non-linear variations of the sky during each scan. To further 
improve the sky-noise removal we hav e used the SURF rou- 
tine calcsky jjenness & Lightfootll998l) . which is similar to the 
remsky routine used for non-raster data. A detailed discussion 
of these methods is beyond the scope of the paper, but while 
scan_rlb relies on the fact that there is no source signal at the 
edges of the map, calcsky makes use of the fact that (a compo- 
nent of) the sky noise is coiTelated across the bolometer aiTay 
during each integration step. The application of both methods 
improved the flatness of the background without changing the 
overall appearance of the image. 
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Table 1. Galaxy parameters 
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Fig. 2. Top: Fitted exponential disk subtracted from the to- 
tal 850 /im emission as shown in Fig. [2 Contour steps are 
6 mJy/beam (1.15 MJy sr~^). The maximum intensity of 
the disk is 66 mJy/beam with a major axis scale-length of 
1.95' . Bottom: Contours of 850 /im in M 51 after subtrac- 
tion of the expo nential disk superimposed on a 6400A opti- 
cal image (:Tilanu s"et alJl988i) . Contours levels are 12 (=1.5 cr; 
2.3 MJy sr~^), 24, 36, 56 and 72 mJy/beam. The image shows 
the centroid of the 850 /im emission to be located along the 
inner edges of the spiral arms and the optical dust lanes. 



Nevertheless, even with the application of the above, rela- 
tively objective, methods background artifacts remained, likely 
because of the missing Fourier zero. This is a common fea- 
ture of SCUBA scan maps. We used an iterative procedure to 
correct the image which involved subtracting an incUned expo- 
nential (the diffuse disk), un-sharp masking and blanking of the 
remaining source signal (arms, nucleus), followed by a polyno- 
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Notes t o Table 1: 

" RSA jSandaee & Tammanrll987h 
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Hc Hl994 



"'Turner & 

" TuUv 1 197 ^ 

Sandaae & Tammannl h975h 

" Beckman et al..t.l99( ) 

^ Laurikainen & SaloH200lh 



mial plane fit and convolution to obtain a smooth, large-scale 
representation of the background. The operations were done 
using the finished coa dded image. While comparable to the 
procedure used by e.g. lAlton et al ] (|2002'), it is subjective in 
nature. As guidelines we used the following criteria: the flat- 
ness of the background and the assumption that the inter-arm 
emission originates from the diffuse disk only. 

The exponential diffuse disk will be discussed separately in 
section 4: we fitted a scale-length of 1 .95' , perha ps fortuitously 
the same as found for the 20-cm 'base' disk ( Tilanus e t alJ 
IT988.) . the inclination and position angle as in Table ^ and 
fitted a peak of 66 mJy/beam. Across the disk of the galaxy 
the fitted background approximates a linear ramp of —12 ± 
10 mJy/beam: from —2 in the south to —22 in the north. The 
mean level of — 12 mJy/beam corresponds to (minus) 9% of the 
132 mJy/beam peak in the raw map. By inspecting the resulting 
inter-arm regions (Fig.|3 we estimate the local background to 
be accurate to a level of about ±2 mJy/beam or about ±5% of 
the typical knots in the arms. Given the overall rms in the image 
of 9 mJy/beam, this uncertainty is about four times larger than 
the statistical uncertainty of the mean over the areas inspected. 
Fig lDshows the calibrated observations after correcting for the 
low-level large-scale background fluctuations. 
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2.2. Total flux density 

Using an aperture of 8' , which just excludes the nucleus of 
NGC5195, we find an integrated flux of 14.5 ± 1.1 Jy. The un- 
certainty quoted indicates the range of values found when vary- 
ing the aperture from 7.4' to 10' , not the uncertainty resulting 
from the background variations. To estimate the latter we in- 
spected the variation of the mean level in 2.5' regions around 
M 51, similar to the size of the nucleus and inner arms. This 
variation is about 0.5 mJy/beam and can serve as an indication 
of the quality of the background correction where it applies 
to regions already largely free from emission. By contrast we 
note that the ±2 mJy/beam uncertainty of the background in 
the inter-arm regions, as mentioned in the previous section, is 
positive in the south and negative in the north, suggesting that 
the uncertainty of the zerolevel taken over the whole of the disk 
of M 51 is less than 2 mJy/beam. Based on these considera- 
tions we adopt an overall uncertainty in the background level 
of 1 .5 mJy/beam corresponding to 1 Jy when summed over the 
8' aperture. Finally we need to add a 10% uncertainty in the 
flux of the calibrators at submm wavelengths. Treating the var- 
ious uncertainties as statistical errors, we derive a total 850 /im 
flux for M 51 of 15 ± 2 Jy. This value includes the contribution 
from line emission within the passband of SCUBA's 850 /im 
filter, a topic which will be discussed in the next section. 

3. Contamination by J=3-2 CO line emission 

The broad passband of the SCUBA instrument includes the 
wavelength of the J=3-2 ^^CO line. Thus, at any point in the 
map the continuum emission measured by SCUBA is, in princi- 
ple, contaminated by CO line emission. In order to investigate 
the extent of the contamination we obtained with the JCMT in 
spring 2001 a set of CO(3-2) spectra across both inner arms 
of M 5 1 and supplemented these with archived spectra of the 
nuclear region. In the arms, the CO line contribution can be 
up to 30-40 percent of the observed continuum emission after 
subtraction of the exponential disk. Thus, a substantial frac- 
tion of the arm morphology seen in the SCUBA map directly 
results from contaminating line emission. Because of the ap- 
parently close corr espondence between CO and dust emission 
(iGuelin etalll 1995 1. the effect is mostly quantitative and does 
not change the overall morphology of the arms qualitatively. 
The J=3-2 emission from the inter-arm regions is weak and 
in individual pointings does not exceed the noise in the mea- 
surements. In the center of M 51, we find an upper limit of 12 
per cent (3(t) to any line contribution to the emission from the 
fitted disk. Consequently, we may assume either a constant CO 
contribution at any level below 12 per cent, or a maximum con- 
tribution of 12 per cent at the center decreasing to zero at the 
outer edge. These two possibilities thus define lower and upper 
limits, respectively, to the 850/xm emission scale-length. 

4. The diffuse disl< 

We have separated emission in the diffuse disk from arm emis- 
sion by subtracting an inclined exponential disk and making 
the assumption that the inter-arm emission originates from the 



diffuse disk only (see also section 2.1). The best fit to the data 
is given by an exponential disk with a maximum intensity of 
66 mJy/beam, and a major axis intensity scale-length of 1.95' , 
corresponding to a linear scale-length of 5.45 kpc at the as- 
sumed distance of M 51. As 1 MJy sr^^ = 5.22 mJy/beam, the 
maximum intensity 66 mJy/beam = 12.7 MJy sr^^. The sub- 
tracted exponential disk and the arm emission are shown in 
Fig. 12] The global parameters adopted for M 51, such as the 
distance, inclination, and position angle, are given in Tabled 

The apparent existence of a diffuse exponential disk raises 
further interesting questions. What kind of dust distribution 
can account for such a disk: an exponential dust distribution 
at a constant temperature, a constant dust distribution with a 
temperature gradient, or both an exponential distribution and a 
temperature gradient? What is the dust mass in proportion to 
the amount of gas? 

4. 1. The Monte Carlo radiation transfer code 
RADMC 

To answer these questions, we have modeled the observed 
emission using a modified version of the axi-symmetric Monte 
Carlo code RADMC (see DuUemond & Dominik i2004) . This 
continuum radiative transf er co de is based o n the method of 
iBiorkman & WoodI tOOlh and ILucvI ( Il999l) . but adapted to 
suit the purpose of this paper. The basic idea of this method 
is as follows. The input luminosity of the system comes from 
the population of stars in the galaxy. At every location in the 
galaxy this is represented by a wavelength-dependent source 
function, with the spectral shape of the local stellar popula- 
tion. Integrated over wavelength and over space, this source 
function yields the total intrinsic bolometric luminosity of the 
galaxy. Each spatial cell represents a certain fraction of this lu- 
minosity, dependent on the spatial distribution of the stars. And 
within such a cell each wavelength interval represents again a 
fraction of this, dependent on the local spectral shape of the 
stellar population. In this way each space-wavelength cell rep- 
resents a well-defined fraction of the total input energy into the 
problem. 

The total input luminosity is now evenly divided over A^phot 
photon packages. These packages are launched into the galaxy 
at random angles from randomly chosen space-wavelength 
cells, where the probability for each space-wavelength cell to 
launch this photon is proportional to the fraction of the total 
input luminosity represented by that cell. Each photon package 
is launched and performs a random walk through the computa- 
tional domain until it eventually escapes to infinity. After that, 
the next photon package is launched. 

As these photon packages travel through the computational 
domain they can experience scattering and absorption events. 
A scattering event merely changes the direction of the photon 
package, but an absorption-reemission event can also change 
its wavelength. The probability function for this random wave- 
length is dP{X)/d\ - dBx{T)/dT, where T is the local dust 
temperature at the moment the photon package enters the cell. 
The dust temperature changes during the simulation, increas- 
ing with the number of photon packages that have entered the 
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cell. Each photon package travelling through a cell increases 
the 'energy' of the cell proportional to the energy of the pho- 
ton package, the length I of the path through the cell and the 
opacity of the cell at that wavelength. The temperature of the 
cell is increased after each passage in such a way that the to- 
tal dust emission of the cell Air pdustii\B\{T)dX x Vcdi 
(where Vceii is the cell volume) equals the sum of all contribu- 
tions pdustK\lLtot / Nphot of all photons that have so far passed 
through the cell. 

The increase of the energy of the cell happens each time 
a photon travels through the cell, while a discrete scatter- 
ing or absorption-reemission event (which changes the di- 
rection and/or wavelength of a photon package) occurs only 
at random locations along the path of the photon package. 
This random location is chosen to be at an optical depth of 
T = I log(ran(isccd))| away from the last discrete event, where 
ran(iseed) is a random number between and 1. 

Once all the photons have been launched and have left the 
system, the dust temperature and scattering source function 
at every location in the galaxy has been determined. Modulo 
random fluctuations this should be the correct temperature 
and scattering source term distribution. Now we use a ray- 
tracing code (part of the radi ative transfer package RAD I CAL 
jPullemond & Turollal '2000')) to produce the desired images 
and/or spectra. But to verify whether this temperature and scat- 
tering source term distribution is indeed the correct solution, 
RADMC has been tested in two ways. First of all, it has been ver- 
ified that flux is conserved to a high level of precision (< 1%) 
by making spectra at many inclination angles and integrating 
over wavelength and inclination to obtain the total output lumi- 
nosity. Since the spectra are produced using a ray-tracing code, 
and not by counting the outcoming photon packages (as is the 
case for the Bjorkman & Wood code), this energy conservation 
is not numerically guaranteed, and therefore represents a good 
check on the self-consistency. The code has also been tested 
against many other codes in a 2-D radiative transfer compari- 
son project (Pascucci et al 2004J . 

4.2. Assumed disk structure 

We assumed that the stellar distribution consists of two com- 
ponents, the exponential disk and the bulge. For the bulge we 
used a Hubble profile (.Revnold&.1913: .Hubble. 19.30,) . which 
models the stellar emissivity as: 

C{R,z) = £disk exp ( - — ] 

\ hs Zs J 

+ i:buigc(l + S2)-3/2 (1) 
in which B is given by the expression: 

Re 

where R en z are cylindrical coordinates, -Cbuigc and £disk 
bulge and disk emissivities at the center of the galaxy respec- 
tively, hs and Zg the respective steUar scale-lengths and scale- 



heights. Re the effective bulge radius and b/a the minor/major 
axial ratio. 

We explored both exponential and constant dust distribu- 
tions. The parameterization to describe the exponential dust 
distribution is the same as used for NGC 891 bv Xilouris et aL 

= expf-^-M") (3) 
\ hd Zd J 

with pc being the central dust density, hd the dust distribution 
scale-length and Zd the dust distribution scale-height. A con- 
stant dust distribution was obtained from the same expression 
by putting the scale-length at infinity. 

Extinc tion coefficients were obtained from the model de- 
scribed bv lDraine&T~^ ill 9841) . We considered four different 
silicate-to-graphite ratios in the grains, 0.5:0.5, 0.4:0.6, 0.3:0.7 
and 0.2:0.8. The particle sizes a were set between 0.005 and 
0.25 pm with a size distribution proportional to a^^'^. Our 
modeling did not include either PAH's or very small dust parti- 
cles. Although these particles dominate the mid-infrared emis- 
sion at e.g. 10-20 pm, their contribution to submillimeter emis- 
sion at 850 /xm is expected to be negligible (see e.g. Li & 
Draine l200ll) . 

The spectral energy distribution of the stars was take n 
from the models calculated by Bru zual & Charloj ill993h . 
Specifically, we used the model for a continuously star-forming 
galaxy with a star-forming rate of 10^^° Mq/ji per unit solar 
mass, assuming a galaxy age of 10 Gyr. We verified that de- 
creases in the assumed age to 5 Gyr or even 1 Gyr do not cause 
significant changes in the parameters needed to fit the 850 pm 
radial profile. The dust temperature rises marginally, causing a 
small change in the dust surface density. The latter, however, is 
much more influenced by the uncertainty in the total luminos- 
ity. 

5. Results and analysis 

5.1. M 51 stellar parameters 

First, the scale-length of the stellar emission must be deter- 
mined as accurately as possible, because it dominates deter- 
minations of the scale-length of the 850 pm dust emission. 
The stellar emi s sivity of M 5 1 has been extensively studied by 
iBeckman et alJ (Il99^ who determined both arm and inter-arm 
scale-lengths in the B, R, I and K photometric bands, where 
they did not correct for extinction. Their inter-arm results are 
summarized in Table[0 Emission in the R and / bands pre- 
dominantly traces old stars which are distributed more or less 
homogeneously through the disk. We used these bands to deter- 
mine the scale-length of the stellar emission. As this emission 
is susceptible to dust attenuation, we have taken extinction into 
account as well. Starting with an initial estimate of the intrinsic 
scale-length of the stars, the radiation transfer code computes 
a value for the attenuated scale-length of the stars, which was 
then compared with the observations. 

The M 51 bulge effective radius was not determined. 
However, as it does not dominate the radial profile, we 
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Table 2. Model results 
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tised a value of 1.03 kpc characteristic for similar galaxies 
llLaurikainen & Salol200ll) . also summarized in Tabled 

Unlike the scale-length of the exponential star distribu- 
tion, the scale-heights of both the stellar disk and the dust 
are in effect degenerate: model and data do not provide useful 
constraints. For instance, a wide range of stellar scale-heights 
yields almost identical 850 fim emission profiles. Dust temper- 
atures are marginally affected, requiring some changes in the 
10 - 100 fim wavelength range of the spectrum. We kept the 
ratio identical to that found in NGC 891 (Xilouris et al] 

11998 V Obviously, as long as the face-on optical depth is kept 
constant the submillimeter emission does not change percep- 
tibly either when we vary the dust scale-height. When we in- 
crease dust scale-heights from 0.25 to 0.35 kpc, temperatures 
change by no more than 1-2 degrees. 

In order to determine the total luminosity of M 51, we 
integrated the UV (OAO), visual and far-infrared (IRAS) 
flux-densities as provided in the NASA/IPAC Extragalactic 
Database, the 170 fim point as given by Tuffs & Gabriel ( 2003) 
and our 850 /im point. We find a luminosity of Ltot = 1.1 x 
IO^^Lq, which corresponds to a starformation rate of a few 
M0yr~^, with an estimated uncertainty of about 30%. Of this 
total luminosity, about 5% originates in the bulge and the re- 
maining 95% in the exponential disk. 
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Fig. 3. The model exponential disk (crosses) compared to the 
actual emission of the M 5 1 disk (solid line), for hd = 7.65 kpc, 
hs = 3.15 kpc, L = 1.1 • 10" Lq, tvj = 4.35 and no CO 
contamination. 
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5.2. Disk parameter determination 

Even when the stellar luminosity and the CO contamination 
would be known, the stellar scale-height, Zs, and the dust scale- 
height, Zd, are degenerate. Adopting values as outlined in the 
previous section, we explored the parameter space over three 
variables: (i) the total amount of dust, which scales with the 
central face-on optical depth ry / integrated from infinity to 
the galaxy center, (ii) the scalelength of the dust distribution 
hd, and (iii) the scale-length of the stellar disk hs, since the ob- 
served scale-length is attenuated by dust. These variables are 
uniquely determined when the stellar luminosity and the CO 
contamination are known. Since this is not the case, we identi- 
fied combinations of parameters best reproducing the observed 
850 fim radial profile. 

We summarized the model results in table 121 An example 
fit is shown in Fig.|3] Note that the model provides a poor fit in 
the very center This is caused by the presence of the bulge in 
the model which provides additional dust heating, hence more 
submillimeter emission in the center We only subtracted the 
disk component from the observed emission. 



Fig. 4. Resulting optical depth for different luminosities as- 
suming either zero CO contamination (upper line) or an om- 
nipresent maximum 12 per cent contamination (lower line). 
The shaded area shows the possible solutions. 

In Figure 13 we show the required face-on central optical 
depth for a range of total galaxy luminosities and CO contri- 
butions. For a luminosity Ltot = 1-1 x lO^^i© these optical 
depths TVJ range from 3.7 (assuming 12% CO contamination) 
to 4.4 (zero contamination). Inclusion of the luminosity uncer- 
tainty increases this to a slightly larger range ry / 3.4^.8. The 
high optical depth derived here is consistent with the typically 
higher values derived from far-infrared and submm c ompared 
to optical observations (see e.g. Bianchi 2004). Popescu et alJ 
( 2000 ) found a value of 3 . 1 for NGC 89 1 . In this case, however, 
arm and interarm emission are fitted simultaneously. The asso- 
ciated scale-lengths vary likewise. For instance, as long as the 
CO contamination is constant with radius, we find a dust scale- 
length hd = 7.65 kpc and a stellar scale-length hs = 3.15 kpc. 
On the other hand, if the CO contaminations is maximal (12%) 
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Fig. 5. Radial temperature profile for an M 5 1 total luminosity 
L=1.11O^^L0. The solid line marks zero CO contamination, 
the dotted line a maximum 12 per cent contamination, and the 
dashed line a contamination decreasing from maximum in the 
center to zero at the M 5 1 outer edge. 



Fig. 6. The dust surface density profiles for a graphite and sili- 
cate ratio of 0.5:0.5 and kappassa — 1.2 g^^cm^; no CO con- 
tribution (solid line), 12 percent CO contribution (dotted line) 
and decreasing CO contribution (dashed line). 



at the center and minimal (zero) at a radius of 10 kpc, we obtain 
a slightly greater scale-length hd = 8.65 kpc. The scale-length 
hs remains the same within the uncertainties.The value of hd 
is somewhat smaller than the scale-length of atomic hydrogen. 
The value of hd is comparable to a scale-lengt h of 9 ib 1 kpc we 
derive for the atomic hydrogen outside 6 kpc jTilanus & AllenI 
Il99ll) . 

In addition to a radial density gradient, we also find a radial 
temperature gradient because the stellar scale-length is much 
less than the dust scale-length, i.e. the stars are more centrally 
concentrated than the dust. Radial temperature profiles, for the 
average grain size at position z = kpc, are shown in Fig.|5] 
When the dust scale-length is larger, we also see a steeper tem- 
perature gradient through the galaxy. The small fluctuations in 
the Fig.|5]profile are artifacts due to statistical noise. 

The good fit of the 850 fim profile provided by exponen- 
tial distributions with finite scale-lengths already suggests that 
the M 5 1 disk is hard to explain in terms of a radially constant 
density distribution. We nevertheless attempted to force such 
a fit, and found that this could only be done if the M 5 1 stel- 
lar scale-lengths were more than an order of magnitude less 
than actually found. Thus, a constant density distribution may 
be ruled out with confidence. An exponential density distribu- 
tion of the cold dust is required in order to explain to observed 
exponential 850 fim disk. 

SCUBA submm observations of the edge-on spiral 
NGC 891 suggest that the cold dust distribution traces the to- 
tal hydrogen column density i.e. molecular hydrogen in the 
central region s and neutral hydrogen at larger radii (Fig. 5 
in Alton et al. 120001). S ubsequent far-infrared observations by 
IPopescu & Tuffs'feOO?) have shown that the cold dust emission 
traced the neutral hydrogen into the extended HI disk. The cur- 
rent observations of M 5 1 , showing an exponential dust density 
distribution, a close association with the molecular hydrogen in 
the central region, and comparable scale-lengths of the dust and 



HI distributions at larger radii further reinforce the suggestion 
that the cold dust is correlated with the total hydrogen column 
density in galaxies. 

5.3. Dust mass and gas-to -dust ratio 

Once the density and temperature profiles are determined we 
may, in principle, determine the radial dust mass distribution as 
well. In Fig.|6]we show the surface mass densities, as a func- 
tion of radius for a graphite and silicate ratio of 0.5:0.5. The 
peak dust surface densities are 0.90 (no CO contribution) and 
0.77 Mq/pc^ (12 percent CO). The central gas (HI + ¥12 + He) 
surface density is Sgas = 60 Mopc^^ (calculated in a future 
paper), with an estimated uncertainty of 25-30%. This results 
in a gas-to-dust ratio of 73±17. 

In this dust mixture, we have a dust emission coefficient 
of K850 = 0.9 g ~^cm~^. Ho wever, the actual value of Kgso is 
rather uncertain. Llames et alJ li2002) used SCUBA observations 
to conclude to a value ksso = 0.7 ± 0.2 g^^cm'^ but their data 
show a spr ead in values for in dividual galaxies by a factor of 
two. Worse. lAlton et alJ ( l2000l) fou nd a three time s highe r value 
for the galaxy NGC 891, whereas 'Agladze et alJ ( 1 199^ in the 
labora tory came up with a value about two times higher than 
that of ljames et alJ ( l2002l) . 

We can make an estimate of ksso in M 51. When we in- 
crease the fraction of graphite in the grains, a higher value for 
the dust emission coefficient is obtained. With a higher dust 
emission coefficient, the surface mass density will be lower. 
For K850 = 1.2 g^^cm^, the peak surface densities as shown in 
Fig.|6]will be 0.65 (no CO contribution) and 0.56 (12 percent 
CO). In table |3] gas-to-dust ratios are shown for four differ- 
ent gas mixtures and for a luminosity of L = 1.1 • IO^^Lq. A 
canonical gas-to-dust ratio of 100±26 is obtained with a dust 
emission coefficient ksr q = 1.2 g ~^cm^, somewhat less than 
twice the value found bv iJames et aLt200Z) . 
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Table 3. Gas to dust-ratios 



silicatergraphite 


K860 (g ^cm^) 


dust: gas 


0.5:0.5 


0.9 


73±17 


0.4:0.6 


1.0 


82±21 


0.3:0.7 


1.1 


91 ±24 


0.2:0.8 


1.2 


100±26 



6. Summary 

1. SCUBA 850 observations of M 51 show a well- 
delineated spiral structure superimposed on a prominent ex- 
tended exponential base disk. 

2. The observed base disk is evidence of an underlying ex- 
ponential density distribution of cold dust in M 51. This 
reinforces the suggestion that the cold dust in spiral galax- 
ies traces the total hydrogen density, i.e. the sum of H2 and 
HI 

3. While throughout M 51 a radial temperature gradient oc- 
curs, because the stars are more centrally concentrated than 
the dust, this gradient by itself can not be the origin of the 
exponential nature of the observed 850 /zm disk. 

4. A reasonable estimate of the dust emission coefficient, 
based on a canonical gas to dust ratio of 100 ± 26 , is 
K850 = 1-2 g~^cm2. 

5. We find a stellar scale-length hs of 3.15 kpc and a dust 
scale-length hd that ranges from 7.65 to 8.65 kpc. 

6. We find a typical central face-on optical depth -ryj = 4.0. 
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